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This work presents a complete cyclic cosmological scenario based on nonlinear magnetic field. 
It is constructed a model composed by five fluids namely baryonic matter, dark matter, radiation, 
neutrinos and a cosmological magnetic field. The first four fluids are treated in the standard way and 
the fifth fluid, the magnetic field, is described by a nonlinear electrodynamics. The free parameters 
are fitted by observational data (SNIa, CMB, extragalactic magnetic fields, etc) and by simple 
theoretical considerations. As result arises a cyclic cosmological model which preserves the main 
successes of standard big bang model and solve some other problems like the initial singularity, the 
present acceleration and the Big Rip. 



I. INTRODUCTION 

Although the inflationary scenario still maintains its 
status of the paradigm in Cosmology, there is an in- 
creasing interest in the alternative proposal of models 
displaying a bounce [THS]. This means the possible ex- 
istence of a collapsing era prior to the actual expanding 
phase. One of the main differences of these two proposals 
concerns the behavior of small perturbations and its evo- 
lution from a perfect spatially homogeneous phase to the 
production of inhomogeneities. There are hopes, due to 
a through analysis of these processes [5], that this ques- 
tion will be settled in the near future. Another important 
point in favor of bouncing models is that they avoid the 
problem of initial singularity. The present paper exam- 
ines one of these scenarios which the dynamical history 
of the universe is partial controlled by a magnetic field. 
It is analyzed its main effects and its compatibility with 
actual observations. The main hypothesis of this model 
concerns the non-linearity character of magnetic field. 

The effects of a nonlinear electromagnetic theory in a 
cosmological setting have been studied in several articles 
[7HH] and for different reasons. For example, in [7] the 
non-linearity is responsible to avoid the initial singularity, 
and in [8 the non-linearity is accountable for generate the 
recent acceleration. In these two papers the framework 
is the same: a cosmological magnetic field governed by 
a nonlinear electromagnetic theory is responsible by the 
desirable effects. 

Recently, these two effects were combined generating 
a cyclic cosmological toy model [TS] . The purpose of the 
present paper is improve this model and investigate how 
realistic it can be. It is a well-known fact that the stan- 
dard cosmological model unavoidably leads to a singular 
behavior of the curvature invariants in what has been 
termed the Big Bang. On the other hand, there are some 
evidences that point in the direction that the universe is 
undergoing an accelerated expansion |161 117) . In princi- 
ple, these two problems have no simple connection and 
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thus does not have a unique combined solution. It should 
be tempted to try to unify these two questions in a sin- 
gle model. It will be shown that in the framework of a 
magnetic universe these two problems are solved at once. 

In general, the cyclic cosmological models can be di- 
vided in two classes: Those which generate the cyclic be- 
havior with non-conventional matter fields [TBI HHj , and 
those which produce the cyclicity through of extension of 
General Relativity 20 22 . The model proposed here be- 
longs to the first class where the non-conventional matter 
is represented by a non-linear magnetic field. 

The plan of the paper is as follows. Section II presents 
the notion of the Magnetic Universe and its generic fea- 
tures. In section III it is fitted or constrained the five 
independent parameters presents in this models. The 
section IV presents the complete scenario consisting of 
the five eras. The paper ends in section V with some 
generic final comments. Moreover, it is presented in the 
appendix A a brief review about the average procedure 
and in appendix B a discussion about the model stability. 



II. MAGNETIC UNIVERSE 

To construct a realistic magnetic universe we start as- 
suming that the universe is composed by the five fluids 
namely baryonic matter, cold dark matter, neutrinos, ra- 
diation and a cosmological magnetic field. The first four 
fluids will be treated in the standard way. Thus bary- 
onic and dark matter are modeling by non-relativistic 
fluids (zero pressure), and neutrinos and radiation are 
featured as ultra-relativistic fluid (The mass of neutrinos 
is neglected). The fifth fluid, the magnetic field, will be 
described by the Lagrangian for the non-linear electro- 
dynamics given by 



Lnled = a F 



/r 
F 
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where the dimensional constants a, (3 and \i are to be 
determined by observation. 

The nonlinear terms can be interpreted as a phe- 
nomenological approach. But a more interesting and 
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fundamental scenario is suppose these terms represent a 
classical interpretation of vacuum polarization. Indeed, 
Heseinberg, Schwinger and others showed that quantum 
corrections due vacuum polarization changes the classical 
Lagrangian introducing nonlinear terms [231 124) . 

Let us suppose that the five fluids are independent. 
Moreover, the magnetic field enter in cosmological sce- 
nario through the average procedure described in ap- 
pendix A. By construction, the average of the electric 
part of the cosmological electromagnetic field must van- 
ishes, i.e. E 2 = 0. And in this case, each term of L^led 
becomes cosmological independent in the sense that each 
one behaves as a non-interacting perfect fluid (for details 
see [H]). 

Treating baryonic matter and dark matter as a single 
fluid and neutrinos and radiation as another single fluid, 
the magnetic universe is completely featured for six in- 
dependent components. In this context, the total energy 
density and the total pressure is given by 



Pt = p m + Pur 
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PB4 = -yPB4- (3) 



The real evolution of the universe depends on the val- 
ues of the six constants presents in pi). However, at least 
in a qualitative manner, it is possible to distinguish five 
distinct cosmological stages just looking the form of the 
terms in These five stages are: 

• The bouncing era: this era is governed by psi 
and it dominates in the earliest universe (before 
nucleosynthesis). In bouncing era the scale factor 
reaches a minimum value. 

• The radiation era: this era is governed by pur with 
a(t) ~ i 1 / 2 . During this epoch occurs the primor- 
dial nucleosynthesis. 

• The matter era: this era is dominated by p m with 
a(t) ~ t 2 / 3, . It is expected that the structure for- 
mation occurs in this epoch. 

• The acceleration era: in this period ps3 govern the 
cosmic evolution. The term p+3p becomes negative 
and the universe begins to accelerate. 



• The re-bouncing era: this era is dominated by Pb4- 
The term p + Zp becomes positive once more and 
starts a new decelerated phase. After a while, the 
a(t) reach a maximum and re-bounces entering in 
a collapsing phase. 



In the section |TVj it will be shown explicitly that the 
universe pass through of all the five stages. 



III. SETTING THE SIX CONSTANTS 

In order to analyze the complete scenario it is neces- 
sary to determine six constants where only five are in- 
dependent 1 . Although five independent constants are a 
elevated number of free parameters, this model describes 
all the evolution for the universe, since the bounce until 
the re-bounce. It keeps untouchable the main successes of 
the standard Big Bang model (primordial nucleosynthe- 
sis, CMB generation, etc) and solve some other problems 
like the initial singularity, the present acceleration and 
the Big Rip. 

The six parameters are divided in two distinct classes: 
the first one, represented by a, Bq and puROi acts on the 
primordial universe; and the second one, represented by 
fx, P and Pmo, affects only the recent universe. Because 
this difference, each class must be fixed by independent 
manner. Let's start with the primordial constants. 



A. Constants for the primordial universe 

In this subsection it will be fixed or constrained the 
three primordial parameters namely a, Bq and puro- 
The most simple to determine is puro, so let 's start with 
it. 

The constant Puro represents the nowadays energy 
density for the all ultra-relativistic particles. It is ba- 
sically composed of radiation p 7 o and neutrinos p V Q. For 
radiation, the precisely CMB measurement [25 a shows 
that, 



SttG 

p 7 o = 4.6 x 10 _5 p c with p c = g = 1.9 x l(T 9 /j" 2 erg 

where h ~ 0.72. 

Theoretical considerations for a massless neutrinos [25] 
imply in 
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1 In this work it is considered only Euclidean spatial section, i.e. 
K = 0. 
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Therefore, the nowadays ultra-relativistic energy den- 
sity is 



Puna = 7.7 x 10- 5 p c 



(4) 



The next parameter it should be fix is Bq. As could 
be seen in ([3]), the constant Bq acts directly only in the 
primordial universe through pbi- Nevertheless, its nu- 
merical value is necessary to determine the a, (3 and p 
constants. Thus, Bq affects not only the primordial uni- 
verse but also, in the indirectly way, the present universe. 

Physically the parameter Bq is related with a cosmo- 
logical magnetic field. Supposing the existence of a large 
scale magnetic field, its background mean value Bq must 
be B < 10 -9 G [27]. In the cgs system, the energy 
density related with this magnetic field is 



= — = 4 x 10 _ ^V era, 



(5) 



where b is just a parametrization (b < 1). 

On the other hand, the present energy density for the 
cosmological magnetic field PB2 is given by 
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So, the comparison between the last equation and ^ 
results in 



where the scale factor was converted to temperature 
through of the relation a ~ T~ l . 

Using T 7 o ~ 2 x 10~ 10 MeV and supposing that the 
nucleosynthesis took place around 1 MeV we finally ob- 
tain 



2 9.5 x 10" 24 
a < p Pc . 



(7) 



For consistence with second condition the bounce 
should happen below the Planck scale, i.e. T^Bou "C 
10 -19 GeV. In this model, the bouncing occurs when 
Pur + Pbx — which implies in 
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and hence 
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At last, using ^ and ([8]), it is obtained the range of 
validity for the parameter a: 
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Replacing ^ and (21) in pi) a first relevant result 
arises: pB2 is always subdominant, i.e. the Maxwell term 
concerning to cosmological magnetic field is not impor- 
tant to the universe evolution. In fact, even the actual 
cosmological experiments including CMB anisotropy ob- 
servation performed by WMap satellite do not have accu- 
racy to measure a cosmological magnetic field. Perhaps, 
the future experiments of CMB anisotropy (e.g. Planck 
satellite) are able to detected it. 

The last parameter necessary to fix is a. This param- 
eter rules the bouncing stage, hence it is concerning at 
the very early universe. These features make difficult 
to precisely determine a. Therefore, instead of to fix it 
we just establish a range of validity for a. The maxi- 
mum and minimum allowed values are obtained through 
of two physical conditions: the bounce must occurs con- 
siderably before primordial nucleosynthesis but it should 
not extrapolate the validity of classical gravitation. The 
first condition implies in \p UR (a nuc )\ < p UR (a nuc ) i.e., 
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Replacing Q and ^ in the last equation it follows 
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The term a 2 ^ 2 at the Lagrangian could be inter- 
preted as an one-loop quantum correction in the infrared 
regime. Indeed, it correction was obtained first by Eulcr- 
Heisenberg [23] • In this context, a 2 values 
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Comparing ( 10 ) with (|9| it could be seen that the both 
results are consistent. So, in this scenario, the bouncing 
stage was generated by quantum electrodynamics effects. 



B. Constants for the recent universe 

In this subsection, we intend to determine the three 
remain constants namely p, (3 and p m o- Differently for 
the three previous parameters, these constants are rele- 
vant at the same time (present universe) . For this reason, 
they should be set simultaneously. 

The procedure that was adopted is based on the fit 
of the model using supernova data. This procedure con- 
cerns only to recent universe, so it could be neglected 
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the terms pbi, Pur and PB2- In this case, the luminosity 
distance g?l is written as 



d L (z) 
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where the Hubble function H(z) is given by 
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with the three Si's defined as 
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(13) 



For 2 = 0, (12 1 implies in O,„o + fl^o + ^/30 = 1- Thus, 
for the present universe, the model has only two free 
parameters. 

The set of SNIa which was used in the fit is the Union 
set [28]. This set collects the mainly sets of supernovas 
which was obtained during the last twelve years. It is 
constituted of 414 SN divided in 13 sub-sets where all 
of them were analyzed through of the same procedure. 
Based on this procedure, the authors excluded 107 SN, so 
only 307 supernovas were regarded true standard candles. 
It is with these 307 SN which is fitted the two parameters. 

The supernova data are usually assumed as Gaussian, 
hence it could be used the statistical likelihood method to 
determine the free parameters. The standard procedure 
consist in to build and minimize a x 2 through the fit of 
the model's parameter. Following the steps point out in 
[2"8] . let's start the construction of x 2 regarding only the 
statistical errors: 
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p,(Zi\fl m O,^fio) = 5 log (4 (Z<|0 m o, O^o) 



- M + a (si - 1) - /3ci 



(15) 



where m., is the B-band peak magnitude, M is the ab- 
solute magnitude, Si represents the stretch in supernova 
lightcurve, Cj represents the color variation between B 
and V bands, a is the stretch parameter, (3 is the color 
parameter and M is a additive factor which include con- 
stants like Ho- The Cjk matrix is the statistical covariant 
matrix, a tot is related with astrophysical dispersion and 
(Ji n t represents the intrinsical dispersion. 

The term Yljf. Cjuej^k represents the uncertainties as- 
sociated with rrii,Si and Cj where ej = (l,a,—/3). In 
principle a and (3 are free parameters which must be fit- 
ted for each specific cosmological model. Nevertheless, 



the authors in 28J found out that a and j3 are almost in- 
dependent of cosmological parameters. So, it is assumed 
that a and j3 are constants 2 , the uncertainties of mi^Si 
and Ci could be directly passed to fi r . . In this case, the 
X 2 is written as 



X 2 (O ro0 ,fi M o,M) =E 



V (p ri - pi (zi\Cl m0 , O^o) - M) 2 



° 2 ot 



(16) 

where the parameter M was embedded in Ai. In 28 is 
available the list with z, p Ti and <r Mr . for each one of the 
307 SNIa. 

Beside a fJir . , there are other three gaussian sources of 
uncertainties: the first one is the dispersion due gravita- 
tional lensing, the second one is the uncertainty in the 
Milk- Way dust extinction correction and the last one is 
the uncertainty due peculiar velocity in the host galaxy. 



• Gravitational lensing decreases the mode of the 
brightness distribution and causes increased disper- 
sion in the Hubble diagram at high red-shift [2"B1 - 
|3"D] . The both strong and weak leasing effects en- 
gender a dispersion of 0.093 x z in magnitude [30j . 
and if the statistical of SNIa is large enough this 
uncertainty could be considered as a Gaussian er- 
ror. 

• All the SNIa light-curve data were corrected for ex- 
tinction caused by our Galaxy. Nevertheless, these 
corrections contain statistical and systematic er- 
rors which act directly on the magnitude value. 
These statistical uncertainties will be neglected be- 
cause three reasons: they only concern the SN with 
z < 0.2, and these ones represent approximately 
20% of total set; out of Galactic plane these uncer- 
tainties value as a rule ~ 0.01; the implementation 
of these uncertainties is very complicated because 
the extinction and its errors change with the Galac- 
tic coordinates. The systematic uncertainties are 
treated later on. 

• As point out in several papers [THl [T71 [551 ISI] > the 

peculiar velocity in the host galaxy causes a error 
a v which varies between 100 and 500 Km/s. In 
this paper, we adopt a v — 300 km/s. So, the un- 
certainty propagated to magnitude due the peculiar 
velocities is 



d h In 10 



2.172 cr v 



where a d = ^ = 0.001. 



2 This approximation engender an new source of error which will 
be regarded as a systematic error. 
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These three sources of uncertainty are encapsulated in 
atot- Thus, the discussion above results in 



(0.093z)' 



2.172 a v 

z c 



(17) 



The last source of statistical error which was consid- 
ered is the uncertainty produced by an intrinsical disper- 
sion <Ji nt ■ This intrinsical dispersion is related about non 
corrected and/or unknown errors. 

To estimate a\ n t it is performed a first fit with a small 
initial value (~ O.Olmag). Then it is determined a new 
Oint imposing \ 2 = 1- And the last step consists in per- 
form the same routine once more to obtain a more accu- 
rate and definitive Oi n t- During this procedure it is not 
used a to t because the intrinsical dispersion is related only 
with the "main" data namely {z,/i,<r M }. 



1. Systematic uncertainties 

The systematic uncertainties are introduced at the 
same way which was performed in |28j . The method con- 
sists in to introduce a gaussian distribution with a AM, 
parameter for each present systematic uncertainty <tam s • 
The new parameters get in additively in the modulus of 
distance, i. e. 

fl (z\Sl mQ , O m0 ) -> ft (z|fimO) + AMs ' 

s 

and the previous likelihood changes to 

Cprev (^mOi Qfj,0) ~ * £post (^mOi ^/jOj AM s ) 

where 



■prev (n m o, fi M o) JJ e 2 " AM * . 

s 

In this kind of construction, we are supposing that a 
specific source of systematic error leads to a "variation" 
in the modulus of distance. And the importance of its 
variation is determined through of a gaussian distribu- 
tion where the standard deviation is identified with the 
systematic uncertainty. 

At the Union set, there are two kinds of systematic 
errors: one which depends on the specific sub-set (e. g. 
due to observational effects) and the other which con- 
cerns to all set of supernovas (e. g. due to astrophysical 
or fundamental calibration effects). The systematic un- 
certainties common to all samples can be absorbed in the 
definition of absolute magnitude. On the other hand, it 
is expected a difference between the nearby (z ~ 0.05) 
and distant (z ~ 0.5) supernovas. So, we can cast the 
common systematic uncertainties to a uncertainty in the 
difference AM between absolute magnitudes of close and 
distant SNs. Following the steps discussed in [23], it is 
defined zjiv — 0.2 as the point which splits the SNs in 



nearby and farther objects. Beside AM, we introduce a 
set of 13 parameters AM S which represents the system- 
atic uncertainties for each one of the 13 sub-sets. 
The complete likelihood for the model is given by 



£post ^ 



(18) 
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The Union set was split in the sum of the 13 sub-set 
whe re eac h su b-set has N s samples. Thus, comparing 

13 

( 16 1 with (19) it follows that J2 N s = N - 



s=l 



The systematic errors arise from seven distinct sources: 

• Errors due the fixation of a and (3 parameters. 

• Errors due the possible samples contamination. 

• Errors due the models of supernovas light curve and 
i^-corrections. 

• Error due photometric zero point calibration. 

• Error due a possible variation of Malmquist bias 
with the red-shift. 

• Systematic errors due gravitational lensing effects. 

• Systematic error due the normalization in the cor- 
rections for Galactic extinction maps. 

Taking account all these effects, the authors of 
reached the following uncertainties: 

(Jam = 0.04 and <tam„ — 0.033 for all s. 

So, the error related with each one of the 13 sub-sets are 
the same. 



2. Results 

First, we take account only the statistical uncertain- 
ties. Using the data presents in [2H], it is performed a 
numerical calculation to determine the cosmological pa- 
rameters. The procedure could be summarized in five 
steps: 

1. Construction of C (£l m o, Q^o, M) not taking ac- 
count a to t- 
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2. Marginalization of A4 parameter 3 . 

3. Calculation of <Tj„ t . 

4. Construction of the new likelihood taking account 
a tot and the a int . 

5. And finally, the minimization of x 2 (^md ^o)- 
The results are shown in figure [T] 
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FIG. 1: Parametric contour plot for Q m o x fi/xo taking account 
only the statistical uncertainties. From inside to outside the 
contour lines correspond respectively to 68%, 95% and 99% 
confidence level. The dashed line represents Q,go = 0. 

The best fit for S! m0 or f2„o could be established 
marginalizing in all the other parameters. Performing 
these marginalization we get in the following results: 



0.44] 



0.03 



and Q 



Ii0 



1.201°;!! • (20) 



Finally, it is determined ttpo through of the equation 
n m0 + f^ + O/30 = 1. Thus, 



do 



-0 fi4~ - 35 
u -°^+0.33 • 



(21) 



The next step is to include the systematic uncertain- 
ties in the procedure. Analogously to the first case, it 
is performed an analytic marginalization at the "non- 
cosmological parameters" namely A4, AM and the set 
of AM S . Again, the method is straightforward but the 
calculation is much more cumbersome than the previous 
case. 

Using the same steps described earlier, it is obtained 
the following result: 



We are interested only in cosmological parameters. 
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FIG. 2: Parametric contour plot for Sl m o x £2/jo taking account 
the statistical and systematic uncertainties. From inside to 
outside the contour lines correspond respectively to 68%, 95% 
and 99% confidence level. The dashed line represents Qpo — 0. 



For this case, the best fits for J7 m o and f2 M o are 

Q m0 = 0.448l°;°39 an d = l-17±o:i? > ( 22 ) 
and the calculated value for fi^o is 

S^^o - -0.62+8:1? . (23) 

As expected, the inclusion of systematic uncertainties 
increases the uncertainties at the cosmological parame- 
ters. 



3. Including a prior in Q, m o 

Measurements of X-ray gas mass fraction f gas done in 
galaxy clusters allow to establish excellent values to the 
ratio flbo/^lmo }32[ I33j. On the other hand, primordial 
nucleosynthesis determine with good accuracy the actual 
barion density flbo- So, join these two facts it is possible 
to obtain an excellent estimation for Q m o- Based on this 
estimation, we intend to build a gaussian prior to fl m o 
as follows: 



P = Exp 



a 



2 < 



where Cl m0 and its uncertainty are obtained through of 
gas mass fraction procedure. 

According to 133) . the f gas is linked with ttbo/QmO as 
shown below: 



f£s DM (*) 



KA^b (z) ( n 



l + s(z) \tt m0 



d\ CDM (z) 
d A (z) 



(24) 
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where d A is the angular diameter distance given by 



d A (z) 



1 



dz 



(l + z)J H(z) 



(25) 



The f^as DM i s the gas mass fraction for KCDM model 
(the reference model) with the following parameters: 
O m o = 0.3, J7a = 0.7 and h = 0.7. The other factors 
which appear in |24|) are described below: 



• K is a constant that parameterizes uncertainties 
in the accuracy of the instrument calibration and 
X-ray modeling. Its value is K = 1 ± 0.1. 

• 7 is related with the non-thermic pressure support 
in the clusters. Its value is 1.0 < 7 < 1.1. 

• b (z) is the depletion factor, i.e. the ratio by 
which the baryon fraction is depleted with re- 
spect to the universal mean. It is modeled by 
b(z) = b (l + a b z) where b Q = 0.83 ± 0.04 and 
-0.1 < a b < 0.1. 

• The function s (z) = so(l + a s z) describes the 
baryonic mass fraction in stars where Sq = 

(0.16 ± 0.05) 4o 2 and -°- 2 < a s < 0-2. 



• The A factor take account the change in the sub- 
tended angle of the clusters due the difference be- 
tween KCDM model and the others models. Ac- 
cording to authors in [33], this effect could be ne- 
glect for all red-shifts of interest, i.e. A = 1. 

To estimate £l m 0i it is used a set of six clusters with 
low red-shifts (see table [I]) 



Cluster classification 


z 


rACDM 1 1.5 
Jgas "70 


Abcll 1795 


0.063 


0.1074(0.0075) 


Abcll 2029 


0.078 


0.1117(0.0042) 


Abell 478 


0.088 


0.1211 (0.0053) 


PKS0745-191 


0.103 


0.1079 (0.0124) 


Abell 1413 


0.143 


0.1082 (0.0058) 


Abell 2204 


0.152 


0.1213 (0.0116) 



TABLE I: The red-shift and f gas for the six galaxy clusters. 
The data were taken from table 3 of reference 1331. 



As we are considering only clusters at low red-shift, it 
can be neglected the dependence on z of the functions 
b(z) and s{z). Thus, inverting (24) it follows: 
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Face on the other uncertainties, the variation of gamma 
could be neglected, i.e. 7 = 1. Using primordial nucle- 
osynthesis models and data of deuterium abundance, it 



is possible to derive the value for Q b o with good accuracy 



n htt h 2 = 0.0214 ± 0.002 with 



0.7 ±0.1, (27) 



where the value of h was determined in |35) . 

The procedure to determine £l m o could be summarized 
in the following steps: 

1. Fit of parameters using the set of SNIa without 
priors. The steps were discussed previously. 

2. Calculation of d A (z) using the previous fit. 

3. Determination of £l m o and its corresponding uncer- 
tainty for each one of the six clusters. 

4. Calculation of weighted mean and the mean uncer- 
tainty for £l m Q. 

5. Finally, it repeats the previous step but now using 
the prior for ft m o. 

These steps was performed some times until ft m o sta- 
bilize in 



il m0 = 0.278 ± 0.038. 



(28) 



The complete result including the statistical and sys- 
tematic uncertainties and the prior in fi m o is shown in 
figure [3j The best fits for f2 m o and fi^o arc 



n m0 = 0.3671°;°^ ari d f^ = 1.47; 



-0.37 

-0.34 ! 



and the calculated value for fi^o is 

^£0 



0.847°; 3 3 I . 



(29) 



(30) 



It is interesting determine explicitly the values of fj, and 
13. Thus, replacing ^ and (|6) in ([l3) follows: 



8.0 x 10 _b 6 p c 



/3 = J-16B^p c n^ = 3.98 x lO^V p 3 J 2 . 

Since p c ~ 10~ 9 erg, the values of p and /3 are ex- 
tremely small and undetectable in the classical earth ex- 
periments. Other interesting feature of this model is that 
the value of fl rn o is greater than the majority models 
presents in literature. Even with the inclusion of the 
prior in f2 m o the difference still remains. It is possible 
that this feature are related with the eminent decelera- 
tion provide by (3 2 /F 2 term. A last point important to 
comment is that with the inclusion of this prior, the re- 
gion where f^o > (below to the dashed line) is excluded 
in 2a. It is a important feature because with positive £lpo 
the model does not engender the re-bounce. 



2.5 
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0.5 




0.3 



0.35 



0.4 



0.45 



FIG. 3: Parametric contour plot for Q m o x fi/xo taking account 
the statistical and systematic uncertainties and a prior in 
fimQ ■ From inside to outside the contour lines correspond re- 
spectively to 68%, 95% and 99% confidence level. The dashed 
line represents f2^o = 0. 



IV. THE COMPLETE SCENARIO 

Before to discuss the complete scenario, it is necessary, 
based on the previous section, choose the values for the 
six constant of this model. For the three constants asso- 
ciated with the primordial universe were selected Q for 



Pc/iWb J|6|) with b — 1 for B 2 and a 



2 = lCT 64 ^ 1 suitable 
with (|9r. For the three constants related with the re- 
cent universe were chosen the best fit for f2 m o, fi^o and 
^/3o given by (29) and (30). This best fit is the most 



convenient choice to characterize the complete scenario 
because it incorporates more observational cosmological 
features than the others two. Nevertheless, using another 
best fit the qualitative results discussed below remain the 
same. 



A. The bouncing period and the primordial 
acceleration 



The conditions for bouncing 

a-Bou = and d Bou > 0, 

implies in 



a Bo 



oO 



OO J QuRO 

which allows rewrite de Friedmann equations as 
O 2 



if, 



2 "r/flo 
o 



OBc 



and 



a 
a 



-H, 



O 2 
o 



(31) 



(32) 



The two previous equations clearly show the existence of 
bounce. Using the numerical values for a 2 , B 2 , and pum, 
it is possible determine the minimum for the scale factor: 



O-Bou 



n 



URO 



1/4 



«0 



7.08 x 10~ 21 



(33) 



where in this calculation the contribution of non- 
relativistic matter was neglected. 



The solution for (31) was studied in [IS] and has the 
following form: 



a(t) = a Bo u\J(2Qtf + 1. 



(34) 



Here the initial condition a(0) = cibou was used. The 
time scale Q is defined as 



Q 



H °S^ =4.1xl0 2 V 1 , 



with H = 72 Km/s.Mpc - see [35]. This solution is 
valid throughout the range in which the non-relativistic 
matter is negligible, e. g. from bounce to primeval nu- 
cleosynthesis. 

The bouncing period plot is shown in figure [4] 



Near to bounce, the Friedmann equations could be ap- 
proximated by, 



and 



where 



V a J V a 



-m 



«0 



uro i — i — 3r2 a0 

a / \ a 



«o 



n Puro , n 16a 2 B^ 
"C/i?o — an<1 "ao — ■ 



Just looking for this plot it is possible to realize 
that the primordial acceleration (inflation) engendered 
by the model is very restricted. In fact, the equation 
( 32 ) states that the acceleration occurs only between 
O'Bou < o,(t) < l.32ciBou, which in usual inflationary lan- 
guage corresponds just to 0.3 e-folds. It rules out the 
standard mechanism to solve de horizon problem. Nev- 
ertheless, in this model the universe is cyclic and eternal, 
so in fact this problem does not exist. Much more in- 
stigate (and complicated) is the question of evolution of 
energy density fluctuations. This issue is directly related 
with model stability (see appendix B) and it will be sub- 
ject of future investigations. 
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FIG. 4: Plot for the scale factor a in function of time t. This 
figure shows that the inflationary phase is very restrained. 



B. The re-bouncing period and the recent 
acceleration 



For the present, the Friedmann equations could be ap- 
proximated by, 



m 



n 



aoV 
' 1 a) 



n 



,,0 



«0 



n 



do 



«0 



(35) 



and 



Hi 



^mO ( — ) 



6n 



v a I 



«o 



ion 



00 



«o 



where Si's are defined in ( 13 ) 



(36) 



The conditions for re-bouncing are given by 
clrbou = and a RBou < 0. 

Unlike the equations of bounce, the equations above 
could not be solved analytically. Even the clrbou could 
not be express from STs in a simple manner. Neverthe- 
less, numerical approaches allow us to study many fea- 
tures of this period. Figure [5] shows the numerical solu- 
tion for (35 1. 



From the numerical calculations involving ( 35 ) it could 



be determined the maximum size for the universe 

aRBou = 1.173a (37) 
and the time remaining to reach it 

tRBou -t = 3.37Gyr. 



The period of acceleration At ace i is obtained from (36) 
and results in 



At a 



5.33Gyr. 



FIG. 5: Plot for the scale factor a/ao in function of time 
t in Giga years, t = is chosen as the instant when the 
actual acceleration begins, to corresponds to interval of time 
between the start of acceleration and the present time. 



It started in 5.31 billions of years ago and in terms of 
red-shift it begun in 

Zacel = 0.533. 

These results show that the acceleration era is almost 
finished, so spending another 20 Million years we are en- 
tering in the re-bouncing period. 



C. The main features for the complete scenario 

To obtain global features for this model it is necessary 
to solve the first Friedmann equation with all terms. It 
is given by 



= m 



BO 



n 



do 



(38) 



where £Ibo — 



is always negligible. 



Again we must use numerical methods to solve ( 38 1 . 
The result is shown in figure [6| 

The plot above shows the complete evolution for one 
and a half cycle. This model is a eternal cosmological 
model, so it is meaningless to talk about the initial time 
and the age of the universe. Nevertheless, one could 
determine the time passed of bounce until nowadays 
Attoday and a period to complete a cycle At cyc i e . Using 
the solution for (38) it follows that At to day = 13.0Gyr 



and A< 



cycle 



32£Gyr 



The solution of (38) allows to establish the values for 
the scale factor at the bouncing and re-bouncing instants. 
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FIG. 6: Plot for the scale factor a/ao in function of time 
t in Giga years, t = is chosen as the instant of bounce. 
The three times ti, t2 and £3 split the expansion period in 
three eras: matter era - 10 -4 < t < t%, acceleration era - 
ti < t < ti and re-bouncing era ti < t < t$. The other two 
eras, radiation and bouncing eras, are too close to origin, thus 
it is not possible distinguish them, at least in theses scales. 



In terms of ao we obtain asou = 7.08 x 10~ 21 a and 
o.rbou — 1.173ao where the maximum variation for the 
scale factor is 1.66 x 10 20 . Comparing these results with 
those obtained previously - equations ( 33 1 and ( 37 ) - it 
verifies that both have the same numerical value. This 



confirms the validity of analysis made in sections IV A 
and lEE 

Another interesting analysis which can be done is the 
evolution of total equation of state lot- Using the relation 



LOT 



Pt 
Pt 



we can study its evolution in terms of time. In most of 
the time during a complete cycle \lot\ < 1- This range 
include two long periods of deceleration (radiation and 
matter eras) and two periods of acceleration (primor- 
dial and recent). However, because of features of bounc- 
ing and re-bouncing, lot must eventually assume values 
less than —1 and values greater than 1. Indeed, near to 
bouncing in a range \t — tsoul < 1-1 x 10~ 21 s we have 
lot < —1, which means a phantom equation of state. 
And near to re-bouncing in a range \t — t^Bou\ < 2.3Gyr 
we have lot > 1, i.e. an Ekpyrotic equation of state. It is 
noteworthy that during all period of recent acceleration 
lot > — 1- This analysis confirm that this model is an 
explicit realization of a quintom cyclic scenario |18[ 136] . 



V. CONCLUSION 

In this paper it was constructed a cyclic cosmolog- 
ical model with five fluids namely baryonic and non- 
baryonic matter, radiation, neutrinos and a cosmic mag- 



netic field described by a nonlinear electrodynamics. 
With five independent parameters it reproduces correctly 
the three expansion phases - radiation, matter and accel- 
erate phases - and produces a bouncing and a re-bouncing 
stages. Another important point is that the main fea- 
tures for the model could be analyzed splitting it in two 
independent phases (primordial and actual) with a link- 
ing phase to connected both. 

In primordial phase the relevant components arc the 
bouncing and ultra-relativistic components characterized 
by the constants a and S1[/.ro- The ultra-relativistic con- 
stant was fixed through CMB measurement and theoret- 
ical consideration about the neutrinos background. On 
the other hand, the constant a remains widely free. It 
happens because there are not much information about 
the pre-nucleosynthesis universe. Another way to de- 
termine a is to interpret the term a 2 F 2 as an one-loop 
quantum correction due vacuum polarization in the in- 
frared regime [53J[23]. In this case, a 2 ~ I0~ 40 p~ 1 and 
this value is in complete agreement with the cosmological 
constraints (see section III A). 

Nowadays the relevant components are the non- 
relativistic matter, the acceleration and the re-bouncing 
terms. The three constants - f2 m oi ^0 an d Qpo - as- 
sociated with each term was fixed using data of SNIa. 
Figures [I] [2] and [3] show that the re-bounce probably oc- 
curs. For the complete fit, represented by figure |3j the 
non-existence of re-bouncing is excluded in 95% confi- 
dence level. Other important feature for this model is 
the high value obtained for fl m o compared with the stan- 
dard ACDM model. Even with a robust prior for non- 
relativistic matter this feature remains. I believe that it 
happens due the eminent deceleration engendered by the 
re-bouncing term. To better clarify this point, it would 
be important to use another cosmological tests. This 
possibility will be investigated in the future. 

A important question concerns about the evolution of 
energy density fluctuations. As point out in section |lV A[ 
the model not engender an inflationary period. Neverthe- 
less, it is not a essential problem for eternal cyclic uni- 
verse since the causal connection could be established in 
previous cycle. Still, the subject about formation and/or 
dissociation of structure (linked with the evolution of 
energy density fluctuations) is an important topic and 
should be investigated in the future. 

Other question concerns about the features of NLED 
used. As mentioned in section [TTJ the NLED was only 
used to describe the cosmological magnetic field while 
the radiation was retracted by the standard electrody- 
namics. So, the non-linear effects act only in the mag- 
netic field and the description of a free radiation remains 
as in Maxwell theory. The main reason to adopt this ap- 
proach is because quantum corrections due vacuum po- 
larization changes the classical Lagrangian introducing 
nonlinear terms |23) . And the new terms only appear 
when we are treating with a quasi-static electromagnetic 
field [21]. So, the fundamental Lagrangian is the Maxwell 
one and the extra terms are just effects of vacuum polar- 
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ization. This approach is directly related with the model 
stability. Indeed, if we consider the Lnled as fundamen- 
tal Lagrangian instead of Maxwell Lagrangian the model 
necessarily will present instabilities near to bounce (see 
appendix B). 

For a wide point of view, we can speculate that the 
electromagnetic phenomenons in extremum conditions 
(e. g. very large and very short scale) are described 
by a non-linear generalization of Maxwell electrodynam- 
ics. Assuming that this generalization in a given scale 
could be expanded in positive and negative powers of 
F = F flv F^ l/ , it is necessary for convergence which 
a k F k > a k+ iF k+1 and (3 k F~ k > f3 k+1 F^ k+1 ^ for k > 0. 
So, it can be argued that in this model the maximum 
and the minimum values for the scale factor limit the 
influence of F k terms with \k\ > 2, i.e. terms of kind 
F 3 ,F 4 ,... and F~ 3 ,F~ 4 , ... are negligible in all cosmo- 
logical time. In this context, the proposal Lagrangian is 
an approximation for a generic non-linear electrodynam- 
ics theory. 
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Appendix A: Review of average procedure 

The standard cosmological model, provided by the 
FLRW geometry, obeys the cosmological principle, which 
means that the 3-space is homogeneous and isotropic ge- 
ometry. Thus, for compatibility with the cosmological 
framework, it is necessary an average procedure in the 
electromagnetic field [37]. As usual, we set the volume- 
tric spatial average of a quantity J by 



J = 



1 

V 



J\J —g d 3 x, 



(39) 



where V = J \J—g d 3 x. And to generate the compatibil- 
ity with the FLRW geometry, it is necessary to impose 
that 



EiEj 



0, B, = 0, EiBi 



1 



0. 



(40) 



1 



B l B J ^--B'g VJ . (41) 

Thus, applying the average procedure in the energy- 
momentum tensor produced by a L = L{F), we obtain a 
T^ v of a perfect fluid where 

?2 



P 



L - - (2B 2 
3 V 



E 2 )L„ 



F ' 



(42) 



and Lp = dL/dF. The bar above F point out that spatial 
average was performed. 



Appendix B: Model stability 

The purpose of this appendix is to discuss the model 
stability. As it is well known all flat cosmological model 
with bounce violates the null energy condition (NEC). 
On the other hand, in [35] was shown that under fairly 
general conditions the NEC violation necessarily im- 
plies instability (ghosts). More specifically, the authors 
showed that when all null vectors violate NEC (as in 
isotropic systems) or when superluminal excitations are 
not present, the NEC violation will result in instabilities. 
These statements extend in a straightforward manner for 
isotropic fluids through the Lagrangian formulation for 
fluid dynamics [39] . Thus, it is reasonable to ask if the 
proposed model is stable? 

To answer this question we will restrict the analysis 
to a region near the bounce. In this case the relevant 
components are represented by the following Lagrangian: 



L 



, o 1 1 

a 2 F 2 ~ -F- -F UR , 



(43) 



where F is associated with the magnetic field and Fjjr 
with the ultrarelativistic constituents. 

To analyze the stability of a model it is necessary to 
perturb this model and analyze the behavior of these 
perturbations. Thus, in the context of magnetic universe 
model, this analysis should be separated in two cases: 
perturbations taken before the achievement of spatial av- 
erages and perturbations taken after this achievement. 

Let's discuss the first case (perturbations taken be- 
fore). The term F 2 arises when quantum corrections 
from vacuum polarization are interpreted classically |231 
[24] . Indeed, Schwinger showed in 24 that for quasi- 
static electromagnetic fields the effects of vacuum polar- 
ization induce corrections in the classical Lagrangian and 
the first of these corrections is given by a 2 F 2 . It is worth 
noting that the complete results obtained in [23] and [24] 
are non-perturbative. So, the term a 2 F 2 represents the 
inclusion (classical) of effects of vacuum polarization in 
the cosmological background. Therefore, the quantum 
perturbations should be done in the usual Lagrangian 
of electrodynamics \F which does not present any in- 
stability problem. In this context, the quantization of 
L ~ a 2 F 2 — jF or any other Lnled is completely mean- 
ingless. 

For the second case, the perturbations must be per- 
formed in energy-momentum tensor of perfect fluid 
(isotropic background). So, the conclusions found in |38] 
can be directly applied, and therefore the violation of 
NEC will imply in an unstable model. Note that these 
perturbations should necessarily be of classical type (such 
as gravitational perturbations) since the spatial average 
has already been taken. In the magnetic universe model, 
this period of instability occurs only near to bounce and 
it quickly ends (see section IV A). 



From a different point of view, we can look at ( 43 ) as a 



fundamental Lagrangian. Thus, it is necessary to analyze 
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the stability of this Lagrangian. The term — \Fjjr will 
be neglected because it does not generate instabilities. 
So, 



2 77.2 



a z F 



1 



-F. 



(44) 



The energy- momentum tensor of Lb is given by 

1 



2 tti2 

a b 



F )g lu , 



remembering that F^ v is determined only by magnetic 
field B. 

To establish if there is a violation of NEC it is necessary 
analyze the quantity T llv n^n v where = (lo, n) is a null 
vector. Performing this calculation in FLRW background 
we obtain the following equation: 



T^n" = [1 - 8a 2 F] 



n x B 



(45) 



where F is function of B 2 . Thus, to happen NEC viola- 
tion (T^n^n" < 0) the first bracket should be negative. 
Another important point is that if T llv n iL n v < then all 
null vector violate NEC. 

In general, the above equation is not spatially isotropic 
since B may depend on points in space. Indeed, it is 



assumed that B has a random spatial distribution whose 
spatial average is determined by (41). Therefore, it is 



possible there are regions of space with instability (where 
NEC is violated for all null vectors) and other regions 
perfectly stable. Note that the instabilities only appear 
at many different points of space if, on average, the term 
8a 2 F is of the order (or larger) of unity, i.e. 8a 2 F > 1, 
and it is exactly what happens near to bounce. So, using 



the bouncing condition (??) and (44) follows 



PBou — -^Fbou — oi 2 F Bou — =>■ 8oi 2 Fbou — 2, 



where F = 2B 2 - eq. (5). 



Since, on average, the signal of ( 45 ) is determined by 



[l - 8a 2 F] or [l - 16a 2 B 2 ] , 

it is possible to conclude that near to bounce most regions 
will present unstable perturbations. However, as the uni- 



verse expands the factor B 



decreases causing a 



decrease in the number of regions with unstable pertur- 
bations. Note that with the increase of a in only one or- 
der of magnitude the probability of having a region with 
unstable perturbations is negligible. One interesting pos- 
sibility for future work is to study how the perturbations 
generated near to bounce evolve. 
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